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We demonstrate two-dimensional robust Raman cooling in a four-level tripod system, in which 
velocity-selective population transfer is achieved by a STIRAP pulse. In contrast to basic 2D Ra- 
man cooling with square envelope pulses [Phys. Rev. A 83, 023407 (2011)], the technique presented 
here allows for a wide variation in the pulse duration and amplitude once the adiabaticity crite- 
rion is satisfied. An efficient population transfer together with attaining of a narrow spread of the 
resonant-velocity group leads to the narrowing of the velocity-distribution spread down to 0.1v Tec , 
corresponding to an effective temperature equal to 0.01T rcc . This robust method opens new possi- 
bilities for cooling of neutral atoms. 



PACS numbers: 37.10.De 

I. INTRODUCTION 

Control of the atomic degrees of freedom at low tem- 
peratures is the starting point for many promising and 
popular research fields which aim either at understand- 
ing and simulating the quantum nature of particle colli- 
sions [1], photoassociation of atoms into molecules [2, 3], 
many-body effects [4] and phase transitions [5] , or at ap- 
plications such as quantum computers [6, 7] and atomic 
clocks [8-10]. Most effects are best observed at ultralow 
temperatures, which can be nowadays achieved for neu- 
tral atoms by a combination of laser cooling and sub- 
sequent evaporative cooling [11]. The latter process re- 
quires also trapping of atoms, usually with a very tight 
confinement, which then precludes efficient use of the 
method e.g. for collimation of slow atomic beams [12- 
14] or in only partially trapped systems such as one- 
dimensional optical lattices [15]. The latter situation is 
utilized e.g. in optical atomic clocks [9, 10]. Although 
tight confinement and the subsequent discrete motional 
state structure for atoms offers many methods for further 
cooling in a manner similar to the cooling of trapped 
ions [16-18] and open possibilities for other interesting 
studies as quantum computing [19, 20] and entangle- 
ment [21-23], alternative approaches are needed in or- 
der to apply cooling at a more general setting such as 
free space. This is the motivation for developing further 
purely light-based methods for reaching similar temper- 
atures as with evaporative cooling, as discussed also in 
our previous work on the topic [24, 25]. 

In the past, powerful cooling techniques have been de- 
signed to achieve subrecoil temperatures of free atoms. 
The "dark state" cooling [26] is very efficient but also lim- 
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ited to rather collisionless situations (low densities). Ra- 
man cooling, on the other hand, is not so density-limited, 
and deep subrecoil cooling in ID has been demon- 
strated [27] and extended to 2D and 3D cooling [28, 29] 
as well. The lowest temperature in 2D, achieved for Cs 
atoms at NIST, Gaithersburg, is 0.15T rcc [29], where T rcc 
is the atomic recoil temperature. The suppression of fur- 
ther cooling is associated with the required cumbersome 
setup of four Raman beam pairs as well as limitations 
of the assumed A-type atomic state system. Our recent 
suggestion of cooling in a tripod atomic level system not 
only reduces the number of Raman beams by a factor of 
two, but also allows one, in principle, to reach temper- 
atures as low as 0.01T rcc . However, more cooling cycles 
are required in 2D Raman cooling as compared with ID, 
which imposes strict demands on the velocity precision 
of the Raman transfer [24]. 

To overcome such a limitation, one can employ the ro- 
bust transfer process provided by STIRAP, as recently 
suggested by us for ID Raman cooling [25]. However, 
the process of transferring atoms collected in the atomic 
"dark state" is not a trivial extension of the ID situa- 
tion, and thereby the 2D case is of special consideration. 
So far, STIRAP in a tripod system by resonant laser 
beams has been experimentally explored only for atomic 
beams [13], although far-off resonant lasers in general are 
used for Raman cooling. This paper demonstrates the- 
oretically 2D Raman cooling by STIRAP going down to 
0.01 Tree, which nevertheless allows a wide variation in 
both the pulse envelope and duration if only the adia- 
baticity criterion is satisfied. The pulse duration needed 
for transfer exceeds the pulse durations for normal Ra- 
man processes, so the advantage of robustness is attained 
only if the cooling time is not a critical factor. This slow- 
ness related to adiabaticity would restrict the application 
of the method in atomic beam collimation to very slow 
beams. Another limitation arises from the specific need 
for a tripod structure, which is not present e.g. at the 
main transitions for the alkaline-earth atoms [30, 31] that 



are currently the strongest candidate for optical atomic 
clocks [9, 10]. 

The organization of this paper is as follows. The nec- 
essary atomic tripod energy state diagram and the corre- 
sponding 2D laser beam configuration are presented and 
discussed in Sec. II. In Sec. Ill we show that, as expected, 
large detuning from the excited atomic state suppresses 
spontaneous decay, and the resonant-velocity group of 
STIRAP under this condition is discussed in Sec. IV. An 
efficient transfer of the atoms from the original ground 
state, even in the case of such a large detuning, leads to 
efficient 2D cooling, for which the parameters are given 
in Sec. V. The cooling itself is investigated in Sec. VI, and 
our research is concluded by the summary and discussion 
given in Sec. VII. 

II. TRIPOD SYSTEM AND LASER 
CONFIGURATION 

Consider a tripod system under conditions following 
closely metastable Ne in Ref. [13]. Pump laser couples 
state (2p 5 3s) 3 Pq to an intermediate state (2p 5 3p) 3 Pi 
(M = 0) which in turn is coupled to magnetic substates 
M = ±1 of (2p 5 3s) 3 P 2 by two Stokes lasers. The pump 
laser is a 7r-polarized running wave propagating along 
axis Oy, and the Stokes lasers are contra-propagating a- 
polarized running waves arranged along the Oz axis (see 
Fig. 1(a)). Note that the metastable Ne system is used 
only as an example. The classical electric field of all three 
laser beams is written as 

E(r,i) = iE P e lkpr - Mpt + e* 5 ^ 1 ^' 

+ ^E s e- tksr -^st + c.c. (1) 

The first term corresponds to pump laser with frequency 
u>p and wave vector kp = kpe y ; the two other terms rel- 
atively correspond to a + - and <j~-polarized Stokes lasers 
with frequencies tjg, uig and wave vectors ks, —kg, 
where kg = kge z . 

Figure 1 (b) illustrates the atomic states coupled by the 
laser configuration, with labelling 

|l) = |(2p 5 3 S ) 3 P ,A/-0), 
|2) = |(2p 5 3p) 3 P 1 ,M = 0), 
|3+> = |(2p 5 3 S ) 3 P 2 ,Af = -l), 
|3-) = |(2p 5 3 S ) 3 P 2 ,A/ = l). 

Taking into account shifts in the centre-of-mass momen- 
tum, we consider an atom of momentum p originally pre- 
pared in state |1, p). Then laser-atom coupling strengths 
are given by 

f|l,p) = ^p(i)e-^ t |2,p + ^k P ), 

V^^ + HkpThks) = ^±(t)e- iw t*| 2) p + ?ikp), 

(3) 
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FIG. 1. (Color online) (a) The 2D laser configuration con- 
sists of three running waves, two circularly-polarized waves 
arranged along the Oz axis and a 7r-polarized wave propagat- 
ing in the Oy direction, (b) The energy-level diagram of the 
tripod system coupled by laser beams, (c) Both Stokes pulses 
fig(t) and Qg(t) forms a counterintuitive sequence with the 
pump pulse Qp(t) in order to carry out atomic population 
from state |1) to |3 + ) and |3~), respectively. 



where V is the coupling operator in rotating wave ap- 
proximation (RWA); the Rabi frequencies 

np(t) = -*^, nf(t) = -^l (4) 

are assumed to be real- valued; d 2 i, djg are the matrix 
elements of the dipole moment operator. Rabi frequen- 
cies (4) evolve in time together with the electric field 
components Ep, E3 and Eg, and have Gaussian en- 
velopes arranged in a counterintuitive sequence as shown 
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in Fig. 1(c): 



n±(t) 



-(t-t P f/2T} 



0± p -(t-t s f/2T* 
S2 s0 e , 



(5) 



where t$ < tp\ here 2Tp, 2T$ are the corresponding pulse 
widths. 

In addition to the atom-field coupling V, the total 
Hamiltonian for the atom-field system 



H = Hn 



2M 



V 



includes the kinetic term P 2 /(2M), and the energy Ho 
of a non-moving atom with the internal state energies 
Ei, E 2 , E + and E$ . As long as spontaneous emission is 
not taken into account, the atomic states in the tripod 
system form a closed family of momentum p: 

F(p)={\l,p),\2,p + Hk P ), 

|3+ p + Kk P - tik 8 ), |3~, p + hkp + hk s )}. (6) 

As a result, in the basis of four bare states, 



ax 
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= exp 
: exp 
■ exp 



Ei 
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2Mh 
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|1,P) 



Ap 1/ 

- A P 



|2,p + ^k P ) 



(7) 



x |3 ± ,p + /ikp T^kc 



the dynamics of the atom is described by the atomic 
Hamiltonian 



H 



/ n P {t) \ 
n P (t) -2A P n+(t) rtg(t) 



nj(t) 25+ 

\ Slg(t) 





2S S J 



(8) 



with the following detunings: 



k P py 
M 



j p , 



M 



(9) 



Here, p^, p z are projections of momentum p on axes Oy 
and Oz, respectively; Ap = wp- W21, A c 



J 23 



the laser detunings; w£ = hkp/(2M), u§ = hk 2 s /(2M) 
are the one-photon recoil frequencies. 



III. SUPPRESSION OF UPPER-STATE DECAY 

The first step of a cooling cycle demands that the con- 
tribution of upper-state decay is as low as possible, be- 
cause the decay broadens the velocity spread of popu- 
lation transfer and thus suppresses the control required 



by subrecoil cooling. To avoid the undesirable sponta- 
neous decay, sufficiently large upper-state detunings Ap 
are commonly utilized. Here such approach means that 
we need to satisfy the conditions 



k PPy k s p z R R 



(10) 



Then the upper state is adiabatically eliminated from the 
Hamiltonian (8) and we can write 



(a 2 |*> 



2A P 



ai|*) 



2A P 



2A P 



oi"l*>, 
(11) 



where I*) is the wave function of an atom. In turn, 
Eq. (11) relies on the adiabaticity constraint 



( a2 |-|*)|«|Ap||( 02 |*)|, 



(12) 



which gives the necessary conditions for the validity of 
the upper-state elimination, namely 



|A P | > Slpo,n%,\6%\ 



(13) 



The latter term shows that the envelopes of the laser 
pulses should evolve in time with a rate that is much 
smaller than the upper-state detuning in frequency units, 
whereas the other terms in the right-hand side respond 
to the splitting of the atomic levels. Then, the reduced 
Hamiltonian in the basis of states {|1), |3 + ), |3 - )} is writ- 
ten as 



H 



2A P 
2J ( 



f2p(i)fi+(t) 



2A P 
V 2h~ P 



26+ 



2Ap 



2A f 



2A P 



2A P 
2A P 



25c 



2A P / 



(14) 



The contribution of spontaneous decay is estimated as 
the loss of population from the upper state of natural 
width r during STIRAP process. With help of the den- 
sity operator a, whose matrix elements are 



,(P) 



1,2,3+3" 



(15) 



the population loss is given by 



d 
dt 



r-Ar 

cr sp (p) = ra 22 (p), cr sp (p) = T / <7 22 (p)d<, (16) 

JO 



occurring during time interval At while the STIRAP 
pulses overlap. Because the adiabatic process takes a 
long time, the overall population loss <7 sp (p) can not be 
neglected at this point. To estimate <7 S p(p), notice that 
the condition in Eq. (11) for upper-state elimination leads 
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to the following inequality: 



4A P 
> .733- (p) <3 



fi+(t) 2 
4A 2 " 3: 



4 3 + (p) 



4A P " 33 



4A 2 P 



(17) 



The effect of the spontaneous decay is now estimated by 

, . n P (t) 2 + n+(t) 2 + n~(t) 2 , . 

Mp) < 3r^^ ^Ar, (18) 



4A p 

and it can be neglected when cr sp (p) <C 1. So, if the 
constraint 



[AH n P {t) 2 + n+{t) 2 + n s (t) 2 

— ^ \Ap~\ 



A7 



(19) 



is satisfied, then the upper-state decay can be neglected 
from consideration. 



IV. ELEMENTARY COOLING CYCLE 

In the first cooling step, STIRAP only accomplishes a 
transfer of atoms through the dark state formed by the 
ground states of the tripod system, thereby determining 
the velocity selectivity of the transfer. As a combina- 
tion of the original ground state |1) with either the |3 + ) 
or |3~) state, the dark state occurs under the condition 
of the two-photon resonance between selected ground 
states. The associated resonant velocities follow from 
the Hamiltonian in Eq. (8) by setting 



5+ = or 



0. 



(20) 



The former condition corresponds to population trans- 
fer by Raman transition |1) ■<-> |3 + ), whereas the latter 
corresponds to the |1) <H> |3~) transition. 

Figures 2(a) and (b) illustrate how the atoms are trans- 
ferred into the |3 + ) and |3~) states depending on their 
velocity, and the corresponding hole burning for atoms in 
state |1) is shown in Fig. 2(c). The use of both transitions 
for transferring the dark-state atoms from the ground 
state |1) is an obvious advantage of the tripod system, 
because atoms can be simultaneously cooled in both di- 
mensions. Form of the burned cross-like hole is defined 
by the laser configuration. The only variable part is the 
position of the cross-like pattern, whose center is given 
by 



4 



0. 



(21) 



and can be shifted by changing detunings Ap, Ag. 

In order to return atoms from states |3 + ) and |3~) to 
the state 1 1} , a fast optical pumping process is needed af- 
ter the STIRAP pulse. The 7r-polarized laser is switched 
off, and only circularly polarized beams are left, being 



(a) 

K4W 




[units of v r , 




[units of v r , 




[units of v r . 



FIG. 2. (a) and (b): Parts of the initial velocity distribu- 
tion transferred by the STIRAP process into the states |3 + ) 
and |3~), respectively; (c) The hole burning arising in state 
|1) due to these transfers. The scaling is given by the re- 
coil velocity v rcc = hkp/M = hks/M and we also set fre- 



^rec 

quency u) R = lo r = uj' s ' 
i± \2 



'* The maximal magnitudes Q, P0 /\A 



p\- 
to 



(ficjn) /|Ap| of the two-photon Rabi frequencies are equal 
2uj r , and A P - A| = u R . The STIRAP-pulse duration T pulsc 
is 96(w H ) _1 . 
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now tuned into resonance. An atom of momentum p is 
excited from either ground state |3 + ) or |3~) to the up- 
per state and therefore it gains a momentum kick along 
the Oz direction. Although the momentum kick equals 
Hu)23/c, one can consider that W23/C ~ ks, and hence the 
atomic momentum becomes p' = p ± hks- Then the 
atom decays into the |1) state emitting a spontaneous 
photon of momentum Ap, where |Ap| = hw2\/c ~ hkp. 
Due to momentum conservation, the atomic momentum 
changes by — Ap, and the population in state |1) takes 
the form 

(l,p|a'|l,p) - (l,p|a|l,p) 

+ (3+ p - hk s + Ap'|«7|3+ p - hk s + Ap') 

+ (3-,p + ftk s + Ap"|ct|3~,p + hk s + Ap"). (22) 

In turn, the populations can be expressed in terms of the 
density matrix elements Cy(p) (15) associated with the 
momentum family J-(p) (6). Taking into account that 
states \dj) (j — l,2,3 + ,3~) are contained in the same 
family .F(p), this leads to the expression 

<r' u (p) = <7 n (p) + (J++ (p - hkp + Ap') 

+ ( 7 33 -(p-^k P + Ap"), (23) 

where momentum shifts are given in relation to ^(p). 

In contrast to the first cooling step where an atom is 
contained in the same family T(p) during the STIRAP, 
the optical pumping process mixes the different families 
J-"(p) as shown in Eq. (23). One can see by averaging 
Eq. (23) over all possible directions of spontaneous de- 
cay that an elementary cooling cycle generally pushes the 
velocity distribution along the Oy axis. If the hole burn- 
ing center v° = p°/M of atoms transferred from state |1) 
(see Fig. 2(c)) is adjusted to u° < 0, then atoms at the 
left-hand wing on axis Oy are pushed closer to the zero 
velocity, which leads to the cooling of the atomic ensem- 
ble. In addition, a laser configuration with a 7r-polarized 
beam in the opposite direction of propagation and adjust- 
ment to Vy > cools atoms also in the right-side wing 
of axis Oy. When these laser configurations alternate, a 
cooling of whole the ensemble becomes feasible. 



V. RESONANT GROUP OF STIRAP 

The efficiency that accompanies the first step of cool- 
ing cycle relies on both a narrow velocity range and the 
entire transfer of resonant-group atoms. In fact, such en- 
tire adiabatic transfer occurs if each laser pulse is tuned 
into resonance with the corresponding atomic transition. 
Such an approach is efficiently applied in Ref. [32] with 
the aim of VSCPT cooling. On the other hand, an effi- 
cient transfer of dark-state atoms takes place even in the 
case of large detuning Ap, as was successfully demon- 
strated for ID subrecoil Raman cooling by STIRAP [25]. 

For velocity-selective STIRAP, a transfer of dark-state 
atoms from the original state |1) evolves with the effi- 
ciency sensitive to the velocity of the dark state. The 



crossing center defined by Eq. (21) is depleted to a greater 
extent, because its velocity v° is attainable for both the 
|1) <-» |3 + ) and the |1) f-> |3~) Raman transitions. For 
the same reason the velocity spread of the resonant group 
is widest at the velocity v°. Next we consider the adia- 
baticity criterion for the population transfer from state 
|1) at the hole burning center, i.e., for conditions given 
in Eq. (21). 

Instead of assuming the conditions in Eq. (21) directly, 
we first take the more general case of Sg = Sg which 
corresponds to an arbitrary velocity projection v y and 



2k s 



Further, we only consider the case of At 



(24) 



when v z — 0. Such an approach gives us a condition 
when the zero- velocity atoms do not leave the |1) state 
and are efficiently accumulated there. 

To simplify the equations of motion, we get from the 
Hamiltonian in Eq. (14) the relationship 



i-(0sW(4|f)-n+(f)(a 3 -|f)) 

= («S+ +i0(t)) (ng(t)<a+|tf) - n+(t)<oi-|*» , (25) 

which only requires that both Stokes pulses, Og(i) and 
flg(t), evolve in time simultaneously: 



"sW = «s(*) = Q(t) 

n+(t) n s (t) 



(26) 



Before the STIRAP pulse starts, the atoms are contained 
instate |1). Hence (a^ | , 5) , (a 3 ~| x I') = 0, and one obtains 
from Eq. (25) that 



ng(t)(at\^}=fl+(t)(a^\^}. 



(27) 



We consider the following coupled (C) and non-coupled 
(NC) states of the coupling operator in Eq. (3): 



|C) = Si,a 



|NC) 



"i(*)|„ + v "sWi 
~W 3; n s (t) 1 



*3 /) 



4 3 h 



(28) 
(29) 



where n s (t) 2 = n+(t) 2 + ng(t) 2 . It follows from Eq. (27) 
that (NCI^) = 0, which in turn leads to relationships 



°^ f) (C|*), (03 j*) = ^4^( c l*>- (30) 



n s (t) 



Equation (26) shows that n^(t)/ris(t) are constant dur- 
ing the STIRAP process. Hence the Hamiltonian in 
Eq. (14) in the basis of states {|ai), |C)} takes the form 



ft (-2S (t) n cS (t) 

2 I Oeff(i) 2(S eS {t) - 6 {t)) / 



(31) 
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The effective detunings and the Rabi frequency are 



VI. FULL COOLING PROCESS 



n P (t) 2 n P (t)n s (t) 

°o{t) = , S2cff(*) = 



4A, 



2A; 



S cS (t) = AS + 



n s (t) 2 -n 2 P (t) 

4A P 



(32) 



where detuning AS determines the offset from the reso- 
nance velocity: 



AS = A s - A P + ^ + w* + u§ = M«v - v° y ). 



(33) 



The Hamiltonian in Eq. (31) describes an effective two- 
level system considered in Ref. [25]. Almost the entire 
transfer of the resonant-velocity atoms occurs once the 
following adiabaticity criterion is fulfilled [25]: 



n P (t) 2 + n s (tf 

\A P \ 



At > 1. 



(34) 



The adiabaticity criterion in combination with Eq. (19) 
gives the conditions 



\a p \ n P (ty 
r 



|Ap| 



•fis(t) s 



At > 1. 



(35) 



This constraint gives the value of At, which is needed for 
achieving an entire transfer of the resonant-group atoms. 
Unlike with the normal Raman processes using square or 
Blackman envelopes for pulses, the requirement of having 
exactly a 7r-pulse is not present here. However, the ap- 
propriate STIRAP pulses should vary slow enough, which 
makes the pulse durations larger than those of 7r-pulses 
in normal Raman process. Hence, the STIRAP transfer 
is suitable in cases where the cooling time is not in any 
critical role, giving in exchange robustness in setting the 
actual pulse durations. 

The resonant-velocity group can be evaluated in terms 
of AS (Eq. (33)). Taking into account that AS = k p (v y — 
Vy) differs from the AS = 2k(v — v ) obtained for the 
ID case [25], the velocity spread of the resonant group 
through v z — is given by [25] 



("£o) 2 + ("g ) 2 
2fcp|A P | 

n 2 

^k P JA p -\ <v ^ v y < 



<Vy-V y < 



O 2 
"po 



2fcp|A P | 



04o) 2 + (»5 ) 2 

2A:p|Ap| 



if A P > 0, 

if A P < 0. 
(36) 



One can see that the velocity spread defined by the two- 
photon Rabi frequencies can get as narrow as needed for 
deep subrecoil cooling. On the other hand, large Rabi 
frequencies broaden the velocity profile of the transfer. 
As a result, an appropriate tuning of the pulse intensi- 
ties allows one to cool the atomic ensemble substantially 
below the recoil limit. 



A single cooling cycle consists of two steps, namely 
the population transfer by STIRAP pulse and the subse- 
quent optical pumping which returns the atoms back to 
the original internal state due to spontaneous decay. We 



assume that k P = k$ and hence ui P 



The 



initial velocity distribution of the atomic ensemble has 
the spread of 3w rec , where v rcc is the common recoil ve- 
locity of all STIRAP lasers. The lasers are detuned from 
the upper state by Ap > 0, and the Rabi frequencies are 
given by 



n 2 P0 = (n+ ) 2 = (n- S0 ) 2 = 2k P \v° y \\A P \, 



(37) 



suppressing the transfer of zero- velocity atoms. The laser 
configuration shown in Fig. 1(a) has the resonant velocity 
of v y < 0, whereas the case of the alternated 7r-polarized 
beam corresponds to v y > 0. After each sequence of five 
cooling cycles with 



\v° y \ = 2 fe -Vcc, 



fc = 0, 



(38) 



the direction of 7r-polarized laser is alternated. The 
broad velocity profiles in the set defined in Eq. (38) in- 
volve all velocity-distributed atoms in the cooling pro- 
cess, whereas those in a narrow velocity group lead to 
the actual deep cooling below the recoil limit. 

The pump and the Stokes pulses have the same pulse 
shape with the pulse half- widths (see Eq. (5)) 

T P = T S =0.5{t P -t s ). 

The duration T pu i se of the STIRAP pulses is defined by 
the start and the end times 



ts — (tp — ts), tend =t P + (tp — ts) 



being equal to 3(t P — ts). As \v y \ decreases, the magni- 
tudes f2p /|Ap| and (flg ) 2 /\A P \ will decrease as well, 
leading to a corresponding increase in At, so that the 
adiabaticity criterion in Eq. (35) is fulfilled. As a result, 
the pulse durations according to the set in Eq. (38) are 
given by 



T, 



pulse 



6 • 2 k r R k = 0, . . 



where t r = (w^) -1 is the recoil time. 

The STIRAP cooling process collects atoms into a nar- 
row peak near the zero velocity, and the height of this 
peak grows simultaneously with the number N of cool- 
ing cycles. The result after applying N = 500 elemen- 
tary cycles of 2D STIRAP cooling is shown in Fig. 3, 
where the height of the peak has become about 230 
times higher than that of the initial broad distribution. 
The velocity spread of the atomic ensemble, given by 
a = (FWHM)/v / 81n2, has been reduced from 3?; rcc to 
0.1v Tec . The corresponding effective temperature T g 
goes down to 0.01T rcc , where T rcc is the recoil-limit tem- 
perature. 
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(b) _ 

Final distribution 



200^ 




FIG. 3. The velocity distribution of an atomic ensemble be- 
fore (a) and after (b) the 2D cooling is applied. The velocity 
spread of (a) 3« rec has been reduced to (b) 0.1v lec after real- 
izing 500 single cooling cycles, corresponding to the effective 
temperature going down to O.OlTrcc- The height of the cen- 
tral peak has increased to about 230 times that of the initial 
distribution. 



Let us consider the result of this cooling as applied 
to metastable Ne atoms under the conditions in the ex- 
periment described in Ref. [13]. Both the <j + - and cr~- 
polarized waves are provided by laser light at wavelength 
As = 588 ran, whereas laser light at \p — 617 nm gener- 
ates the 7r-polarized beam. As a result, the final effective 
temperature of cooled atoms is given by 

T cS « 0m ^^( k P + k s) = 26 nK, 
where fcs is the Boltzmann constant, M is the Ne atomic 



mass. The duration of only STIRAP pulses in our scheme 
is 18600 t r , hence full cooling takes about 0.1 s. 
VII. CONCLUSIONS 

We have considered a variant of optical cooling based 
on velocity-selective STIRAP transfer in a four-level tri- 
pod system. This approach extends into 2D the recently 
proposed ID cooling method [25], providing strong trans- 
verse cooling below the recoil limit. In contrast to the 
normal 2D Raman cooling [24] , the method is robust and 
versatile as long as the adiabaticity criterion is satisfied. 
Strong and efficient cooling is especially attainable at the 
limit of large detuning from the upper state in the tripod 
configuration. The numerical results demonstrate a 2D 
cooling down to 0.01T rec . 

As a topic of discussion, we note that the success of 
evaporative cooling in reaching the atomic phase-space 
density that is required for quantum degeneracy has di- 
minished strongly the original interest in light-based cool- 
ing methods, although the increasing variety of experi- 
ments with neutral atoms under differing circumstances 
is reviving this interest. Similarly, the dynamics of light- 
assisted cold atomic collisions [1, 33] have not been fully 
explored, and in fact the issue of their character is still 
unresolved experimentally [34, 35]. The dynamical view- 
point involving level crossings [36-38] and the comple- 
mentary view of steady-state description [39, 40] have 
both their supporters, and it would be of interest to ex- 
amine the dependence of the collisional atomic kinetic 
energy gain as a function of laser intensity and especially 
detuning [41, 42]. At ultralow temperatures the collisions 
take a very different nature compared to the more semi- 
classical idea of atoms approaching each other [1, 43]; 
also the interesting question about the role of the higher 
partial waves tends to disappear when quantum statistics 
steps in, and energies limit the processes to the s-wave 
only [44]. Such studies are an example where one could 
apply such light-based cooling methods as we have pro- 
posed. A special feature in the STIRAP-based cooling 
is the possibility to use large detunings. This reduces 
the role of light-assisted collisions or reabsorption of scat- 
tered photons in the cooling process, allowing higher den- 
sities than available at standard magneto-optical traps, 
while collisions and other properties of a cold but still 
non-degenerate gas can be analysed with separate probe 
lasers. 
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